Critical Fluctuations and Disorder at the Vortex Liquid to Crystal Transition 
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We present a functional renormalization group (FRG) analysis of a Landau-Ginzburg model of 
type-II superconductors (generalized to n/2 complex fields) in a magnetic field, both for a pure 
system, and in the presence of quenched random impurities. Our analysis is based on a previous 
FRG treatment of the pure case [ E. Brezin et. al., Phys. Rev. B,31, 7124 (1985) ] which is an 
expansion in e = 6 — d. If the coupling functions are restricted to the space of functions with 
non-zero support only at reciprocal lattice vectors corresponding to the Abrikosov lattice, we find a 
stable FRG fixed point in the presence of disorder for 1 < n < 4, identical to that of the disordered 
0(n) model in d — 2 dimensions. The pure system has a stable fixed point only for n > 4 and so 
the physical case (n = 2) is likely to have a first order transition. We speculate that the recent 
experimental findings that disorder removes the apparent first order transition are consistent with 
these calculations. 
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PACS numbers: 74.60. Ge, 64.60.Ak 

In 1957 Abrikosov |J demonstrated the existence of the 
'mixed phase' in type-II superconductors wherein the ex- 
ternal magnetic field penetrates the sample in the form 
of a triangular flux line array which now bears his name 
- the Abrikosov lattice. This work was based on a mean- 
field analysis of a Landau-Ginzburg formulation of the 
problem, and the transition from the normal phase to 
the mixed phase was found to be continuous. The natural 
question which then arises is: how do fluctuations affect 
the nature of this phase transition? This question was 
finally addressed by Brezin, Nelson and Thiaville (BNT) 
in the form of a functional renormalization group (FRG) 
analysis of the Landau-Ginzburg theory Their con- 
clusion, based on a numerical analysis of the (functional) 
flow equation, was that fluctuations drive the transition 
from continuous to first order in spatial dimension d < 6. 

Our original motivation for studying this problem was 
to investigate the effect of quenched random impurities 
on the order of the transition, as it is well known that 
such randomness can often force a first-order transition 
into a continuous one |3j . An interesting example of this 
is the transition from the normal metal to the Meissner 
phase which although first-order for a pure sample j|, 
is driven to a continuous transition in the presence of 
disorder |J. Our analysis was based on the work of BNT. 
In the course of our work, we became aware that the 
representation of the FRG used by BNT might not always 
be sensitive to the existence of fixed points and therefore 
their conclusions on the nature of the transition should 
be re-examined. 

The results from our analysis for the pure system are, 
however, qualitatively the same as those of BNT: no sta- 
ble fixed-point exists in the FRG, which may be consid- 
ered as evidence of a first-order transition. The inclusion 
of quenched random impurities changes this result dra- 



matically. We find that the FRG supports a stable fixed 
point if the coupling functions (to be defined shortly) are 
restricted to the space of functions with non-zero support 
only at reciprocal lattice vectors {G} which correspond 
to the Abrikosov lattice. 

The starting point of the analysis is the Landau- 
Ginzburg free energy functional for type-II superconduc- 
tors in the presence of quenched random impurities 
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where {ipi} are a set of n/2 complex fields (the implic- 
itly repeated indices i and j are to be summed from 1 
to n/2), H is the external magnetic field, and St is a 
field representing the quenched random impurities whose 
quenched average is zero, and whose correlator is given 
by (St(t)St(t')) = A a S d (r - r'). As defined the theory 
has two coupling constants, go and Ao. We have written 
the theory for arbitrary spatial dimensionality d. In the 
physical case of d — 3, the external field H picks out a 
transverse plane (x, y) in which the Abrikosov lattice is 
formed. In general d, we define the field H to be directed 
in a d — 2 hyperplane , so that the Abrikosov lattice is 
still confined to the two-dimensional (x, y) plane trans- 
verse to this field. 

We shall first concentrate on the pure system (Ao = 0). 
Since we are interested only in the critical region, there 
are simplifications which may be made to the form of 
the free energy. We refer the reader to Ref. M for 
the details. The key point is that one may restrict 
ones attention to the lowest Landau level, which cor- 
responds to writing the order parameter in the form 
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"4>i(x, y, rj_) = 4>i(z, rj_) exp(— /i 2 z*z/4) where z = x + iy 
and /x 2 = e*H. A perturbative analysis of this simpli- 
fied theory reveals that one-loop contributions diverge 
due to critical fluctuations in the (d — 2) hyperplane , 
thus leading to an upper critical dimension defined by 
d c — 2 = 4. To control these divergences, one implements 
a renormalization group procedure (perturbatively in the 
parameter e — 6 — d) and, as found by BNT, the theory is 
not closed under renormalization; i.e. the theory defined 
with a simple coupling constant go is not a sufficiently 
general theory. To cure this problem, BNT generalized 
the coupling constant to a coupling function, which un- 
der the influence of translational invariance leads one to 
consider the new bare theory: 

F = J d d - 2 r± _ Jdz dz* e-" S l*l 3 / 2 (|Vx&| a +T|&| a ) 
+ J d d - 2 r ± J dztdzl e~^'^' 2 j dz2<lz * e -M 2 NI 2 /2 
xg(\ Zl -z 2 \ 2 ) |0j(zi,r x )| a Mz 2 ,r ± )\ 2 . (2) 

To close the FRG flow equations, BNT then made a rep- 
resentation of the coupling function g(u) of the form 
g(u) — | J" 00 ^ p(x)e~^~/ 4xS>u leading to a (compli- 
cated) FRG flow equation for the 'weight' p(x). Numer- 
ical integration of this flow equation failed to lead to a 
fixed-point form of p(x) and this was interpreted by BNT 
as signalling a first-order transition. 

One point we wish to make in this Letter is that rep- 
resenting the coupling function in terms of p(x) is rather 
unfortunate. One can show that stable FRG fixed points 
exist for the model described by Eq.(|^) for the cases 
n — ► oo H and n = 0. In each case, the weight func- 
tion corresponding to the fixed-point form of the coupling 
function is ill-defined. In other words, a failure to find 
a fixed point form for the weight p(x) is not necessar- 
ily related to the existence or otherwise of a fixed point 
form for g(u). Details of this rather technical point will 
be given in a future publication (tJ , together with the full 
derivation of the results given below. For now, we shall 
simply indicate how to proceed in an alternative fashion. 
Instead of using the integral representation of g{u) de- 
fined above, one may construct a FRG flow equation for 
the Fourier transform (FT) <?(k) of the coupling function 
(we stress that k is a two-dimensional wave-vector.) It 
then turns out to be more convenient to concentrate on 
the function /(k) = e~~ k / 2ai g(k). By studying the two- 
point correlation function one then finds that the Hartree 
and exchange diagrams (whose sum gives the one-loop 
correction to the correlation exponent) are given respec- 
tively by /(0) and /(0) (where / is the inverse FT of /) 
indicating that this choice of representation is a natural 
one. 

As mentioned above, one can show the existence of 
stable FRG fixed points for the special limits of n — > oo 
and n = 0. We shall briefly describe these two cases. 



In terms of the coupling function / one may derive the 
following FRG flow equation 

dj(k) = ef - (n/2)/ 2 - 2/ o / - 2ff* (3) 

which is in a differential form, with e l corresponding to 
the scale change of the FRG. We have introduced the 
notation: 

a o (3 = J ^ a(p)0(k - p) cos 2 (p x k/2) (4) 

and 

a*(k)^ J a(p)cos(pxk) (5) 

with the definition p x k = p x k y — p y k x (we have also 
scaled wave-vectors by /J, and / by l/(2/i).) 

For n — > oo the flow equation is trivially solved. We 
have dif — ef — f 2 (for the rescaled function / — > 2//n), 
with stable fixed-point solution / = e. This solution in- 
dicates that the coupling 'function' g(u) has the nature 
of a distribution |(| (leading to problems of definition 
for the weight function of BNT.) One may consider 1/n 
corrections to this result by analysing the full flow equa- 
tion, Eq.(^|). One sees that a simple analytic expansion 
in powers of 1/n is not possible, so the limit of n — > oo 
must be seen as somewhat singular. The other extreme 
of ri = may also be considered in a simple way. Setting 
ri = in Eq.(^), indicates that one may obtain the fixed 
point /(k) = (e7r/2)<5 2 (k). A stability analysis shows 
that this fixed point is stable against arbitrary pertur- 
bations. This fixed point corresponds to g{u) = const, 
(which again leads to problems in the definition of the 
weight function of BNT.) It is noteworthy that a naive 
expansion around this result in powers of n is not possible 
- the case of n — is therefore also singular. An inter- 
esting possibility is to expand about the n = limit by 
making a scaling Ansatz of the form / = (l/7i)s(k/n 1 / 2 ), 
where the scaling function s is to be determined from the 
FRG. Such a strategy has so far yielded no useful results, 
and a full discussion will therefore be postponed until a 
future publication M. 

Unfortunately the flow equation for / for the physical 
case of n = 2 is not analytically tractable and one must 
resort to numerical integration. This task is highly non- 
trivial as the flow equation has the form of a non-linear 
integro-differential equation. We have failed to find any 
stable fixed points for this flow equation for finite n > 0, 
although we have been able to generate many unstable 
fixed points by means of a Newton root-finding scheme 
applied to the fixed-point form of the flow equation. Al- 
though we cannot say with certainty that stable fixed 
points do not exist, we are inclined to believe there are 
none. This 'runaway' of the coupling function signals a 
breakdown of perturbation theory which may be inter- 
preted as evidence of a first-order transition. We shall 
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present an alternative scenario for n > 4 later in this 
Letter. 

We now turn to the first main point of this Letter - 
the effect of disorder on the transition. One may proceed 
through all the steps that led to the model described in 
Eq. (Q) , but now with the addition of the disorder field. 
[One uses the usual tricks to average over the disorder 
p|, thus producing an effective quartic interaction with 
strength (-Ao).] In a completely analagous fashion to the 
pure case, one quickly finds that the theory with simple 
coupling constant Ao is not closed under renormalization, 
and one is therefore obliged to generalize the theory by 
replacing Ao with a function A(|zi — Z2\ 2 )- We choose 
to make the natural FT representation as before. We 
therefore have flow equations for the functions /(k) and 
D(k) = e - fe W A(k) of the form 

dif(k) = ef - (n/2)f 2 - 2/ o / - 2/7* + 4/ o D + 2fD* 

(6) 

diD(k) =eD- nDf + 2D o D + 2DD* - 2Df* 

with the same notation as used above. 

We investigated these coupled flow equations numer- 
ically for the physically interesting case of n = 2 and 
found no stable fixed points, but only a runaway of the 
coupling functions. Again, this breakdown of perturba- 
tion theory could be interpreted as signalling a first-order 
transition. In the remainder of this Letter we shall main- 
tain an alternative scenario regarding this breakdown of 
perturbation theory which leads to very interesting con- 
sequences. 

We shall interpret this breakdown of perturbation the- 
ory as a flow to some strong-coupling fixed-point. In 
terms of the model representation in Eq.^, this fixed- 
point is not accessible via an e-expansion. However, we 
may envisage constructing a "</> 4 " model which is per- 
turbatively close to this fixed-point. One would expect 
that any model which is a candidate for describing the 
vortex liquid - vortex crystal phase transition would be 
able to describe the two broken symmetries of the low- 
temperature Abrikosov phase. One broken symmetry is 
ODLRO, or phase coherence in the directions rj_, and 
the other broken symmetry is the breaking of full ro- 
tational and translational symmetry of the liquid phase 
down to those of the triangular Abrikosov lattice. A 
model which can capture these broken symmetries arises 
if one considers the same theory as before but with the 
added constraint of restricting / and D to have non-zero 
weight only at values of k coincident with the set of re- 
ciprocal lattice vectors (RLV) {G} of the Abrikosov lat- 
tice. In practice we set /(k) = ^ G A(G)S 2 (k — G) and 
D ( k ) = Eg B(G)5 2 (k - G). With these forms one sees 
that the new field theory is closed under renormalization; 
coupling functions not of this form are not produced per- 
turbatively. Furthermore, the flow equations given above 
are still valid. Of course such a model cannot describe 



the liquid phase. In that respect it is similar to some of 
the theories of two-dimensional melting which, being ex- 
pressed in terms of Burgers' vectors, can only in a staight- 
forward way describe the crystalline phase (see e.g. Q). 
While we have been unable to derive this RLV model 
starting from Eq.(||), we feel it captures the essence of the 
symmetries broken in the transition and hence because 
of 'universality' arguments provides a way of calculat- 
ing critical exponents etc. at the transition, should it be 
a continuous one. Another noteworthy point is that the 
FRG flow equations in Eqs.@ determine /(|k|), and it is 
hard to see how the rotational symmetry can be broken 
down to that of a triangular lattice unless such break- 
ing is introduced 'by hand', as indeed is done within our 
model. The introduction of the RLV model is the second 
major point of our Letter. 

One finds that the flow equations are immensely sim- 
plified by our particular choice of the RLV - the RLV 
of the triangular Abrikosov lattice. For this choice the 
trigonometric terms in the convolution-type integrals be- 
come equal to unity, and the functions / and D become 
self-reciprocal under the FT. Proceeding with this choice 
of the RLV, it is possible to make analytic progress with 
the simple Ansatz A(G) — A, B(G) = B, i.e. choosing 
the RLV coefficients to be independent of G. 

For this Ansatz one easily finds fixed-point values for 
A and B given by A FP = e/(2(n - l)5 2 (0)) and B FP = 
(4 — n)e/(8(n — 1)<5 2 (0)) which can be proven to be stable 
under arbitrary RLV perturbations for 1 < n < 4. The 
symbol <5 2 (0) is proportional to a free sum over RLV, or 
equivalently to the system size. The one-loop correction 
to the correlation length exponent is given by 2 — 1/v = 
(n/2)/(0) + /(0) + L>(0), and we therefore have at the 
above fixed point 
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To address the question of stability of this fixed point, 
we consider perturbations to the coupling functions. In 
order to remain within the space of coupling functions 
consistent with the RLV model, the perturbations must 
in turn be restricted to the RLV. Explicitly, we write 

/ = f F p + a(G)S 2 (k-G) 

G 

D = D F p + Y b{G)S 2 {k - G). (8) 

G 

Substituting these forms into the flow equations Eqs. (|^) 
and then performing a linear stability analysis yields the 
eigenvalue spectrum for the (infinite dimensional) stabil- 
ity matrix. Stability is ensured by requiring that this 
spectrum has no positive eigenvalues. It turns out that 
the spectrum is characterised by only two different eigen- 
values, a non-degenerate eigenvalue Ai = — e and a highly 
degenerate eigenvalue A2 = e{n— 4)/(4(n— 1)) (to leading 
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order in e), which clearly shows the stability of this fixed 
point in the range 1 < n < 4. Interestingly the value 
of the correlation length exponent and the value of the 
stability eigenvalues (which are related to correction-to- 
scaling exponents) are the same as those obtained for the 
simple 0(n) model in the presence of disorder M (but in 
two dimensions higher; note that in this Letter e = 6 — d.) 

One may also use this RLV Ansatz in the absence of 
disorder, i.e. Eq.(||). In this case one finds a fixed-point 
solution App = 2e/((n + 8)<5 2 (0)) which is found to be 
stable to arbitrary RLV perturbations only for n > 4. 
The value of v at this fixed-point for n > 4 is identical 
with that obtained from the pure 0(n) model Jl(]] (again 
in two higher dimensions.) The stability analysis for the 
pure case reveals an eigenvalue spectrum characterised by 
two different eigenvalues, with values Ai = — e and A2 = 
e (4 — n)/(4(ra + 8)). Clearly the fixed point is only stable 
for n > 4, and one may therefore not make the direct 
connection to the pure 0(n) model, since n — 4 plays 
no special role in that case. (Intriguingly, the value of 
the eigenvalues for the pure superconductor are the same 
as those for the Heisenberg fixed point in the disordered 
0(n) model § .) 

The results for n = 2 in the presence of disorder sug- 
gest that our problem is in the same universality class as 
the disordered 0(n) model, but in two dimensions higher. 
Hence the lower critical dimension of our problem would 
be expected to be 4. That is, in the presence of disorder, 
one would not expect there to be a phase transition in 
d = 3 or d = 2. However, our perturbative treatment of 
the disorder is appropriate only for weak disorder. Strong 
disorder might drive the system to a gauge glass phase. 
However, the recent simulation results of Bokil and Young 
pT[ suggest that the lower critical dimension of the gauge 
glass is greater than 3. The old argument of Larkin |1J] 
suggests at the very least that disorder modifies the na- 
ture of any transition below 4 dimensions. 

In the absence of disorder, we have been unable to find 
for n — 2 any stable fixed points even within the RLV 
model, and deduce that the original conclusion of BNT 
that the transition becomes first order below 6 dimen- 
sions is likely to be correct. Our calculation gives no 
information on the lower critical dimension of the pure 
case. One of us has argued that thermal excitation 
of phase fluctuations does not permit the simultaneous 
existence of ODLRO and the vortex lattice for d < 4. 
However, while the RLV model certainly ensures that 
the low-temperature phase (if it exists) is crystalline, it 
does not necessarily require that phase to have ODLRO. 

Hence, it is very tempting to compare our predictions 
- that in the absence of disorder a first order transition 
should exist and that in the presence of (weak) disorder 
no transition is to be expected - with the experimen- 
tal results of Fendrich et al. [O], who found that they 
could suppress the apparent first order transition seen in 
untwinned single crystals of YBCO by introducing point 



defects via electron irradiation. Futhermore, the resistiv- 
ity in the disordered system did not show the behaviour 
expected according to the gauge glass phase transition 
hypothesis. To our eyes their resistivity curves after ir- 
radiation look consistent with a gradual pinning of the 
vortices as the temperature is lowered, i.e. with the ab- 
sence of a phase transition. 
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